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MAXIMUM  LIKELIHOOD  TRAJECTORY 
DETERMINATION  FROM  RADAR  DATA 


I .  INTRODUCTION 

A  problem  of  considerable  importance  in  the  study  of  radar 
tracking  data  is  how  best  to  obtain  a  target’s  trajectory  taking 
into  account  the  relative  accuracy  of  range  and/or  doppler  meas¬ 
urements  over  the  coarser  angle  measurements.  Various  formula¬ 
tions  of  this  problem  have  been  made  and  this  report  draws  upon 
these  previous  works.  This  report  describes  a  maximum  likelihood 
determination  of  the  unknown  parameters  defining  a  target's  tra¬ 
jectory.  Illustrations  of  the  accuracy  and  speed  of  the  method 
are  made  with  examples  of  fitted  satellite  data  from  the  TRADEX 
radar.  The  procedure  obtains  a  solution  in  3  iterations  or  less. 

Basically,  the  analysis  contains  these  important  elements. 

The  method  is  a  general  maximum  likelihood  method  in  that  it 
obtains  a  reference  trajectory  from  the  data  and  then  perturbs 
the  parameters,  in  this  case,  the  six  initial  conditions,  to 
produce  a  time  history  of  the  resulting  changes  in  the  trajectory. 
These  differences  can  be  viewed  as  a  manifold  of  linearly  inde¬ 
pendent  functions  of  time.  The  data  residuals  are  assumed  to  be 
members  of  this  manifold  and  the  components  which  give  the  proper 
corrections  to  the  initial  conditions  are  uniquely  determined  by 
maximization  of  a  likelihood  expression. 
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Section  II  gives  a  description  of  the  techniques  with  a 
basically  complete  but  simplified  derivation  of  the  maximum  like¬ 
lihood  algorithms.  Section  III  contains  illustrations  of  the 
accuracy  and  fast  convergence  of  the  method  with  curves  and  tables 
of  results  from  TRADEX  radar  data.  An  interesting  study  of  how 
the  trajectory  initial  conditions  deteriorate  as  shorter  data 
spans  are  used  is  given.  Section  IV  contains  some  illustrative 
cases  where,  starting  with  elevation  and  range  data  without  re¬ 
fraction  corrections,  the  method  obtains  residuals  which  roughly 
agree  with  the  normal  values  for  refraction  correction.  Section  V 
presents  the  conclusions  from  this  study  and  suggests  other  areas 
for  further  investigations.  An  appendix  gives  a  derivation  of  a 
refraction  correction  formula  used  in  the  data  reduction. 

II.  DESCRIPTION  OF  TECHNIQUE 

The  technique  will  first  be  described  in  an  abstract  way  and 
then  the  procedure  used  by  the  computational  program  will  be  out¬ 
lined  in  simpler  terms.  In  general,  the  equations  of  motion  of 
the  target  are  supposed  to  be  known,  except  for  a  certain  number 
of  constants  or  parameters  (P^) •  In  the  case  of  the  satellite 
problem,  with  the  gravitational  model  uniquely  defined,  the  only 
parameters  are  the  six  initial  conditions;  and,  henceforth,  we 
set  k  =  1,6.  The  analysis  then  proceeds  along  the  following 
lines : 
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A  reference  trajectory  is  obtained  from  a  set  of  zeroth- 

order  parameters  (P°) •  Any  neighboring  trajectory  has  its 

parameters  incremented  by  a  6-component  vector  (x^),  which  is 

to  be  added  to  (P°) •  The  differences  between  a  neighboring 

trajectory  and  the  reference  one,  in  the  radar  observables  range 

r,  elevation  e,  and  azimuth  a,  are  functions  of  time  t.  Let 

(x.)  ,  n  =  1,6,  be  six  different  increment  vectors  defined  sue- 
k  n 

cessively  by  the  component  arrays:  ( 1 , 0 , 0 , 0, 0 , 0) ,  (0 , 1 , 0 , 0 , 0 , 0) 

. ( 0 , 0 , 0 , 0 , 0 , 1 ) .  One  then  computes  a  set  of  eighteen 

(6  by  3)  "influence  functions",  which  represent  the  perturbations 
6r  (t),  6e  (t),  and  6a  (t),  n  =  1,6,  due  to  the  unit  vector  in- 
crements  in  one  of  the  six  parameters.  The  perturbations  corres¬ 
ponding  to  parameter  increments  given  by  the  array  ( , 

6 

XC,XC)  are  then  given  by  6r(t)  =  E  X  5r  (t),  in  which  the  line- 
5  b  ,  n  n 

n  =  1 

arity  in  the  coefficients  (X^)  is  a  basic  assumption  for  the  maxi¬ 
mum  likelihood  method.  Any  departure  from  linearity  can  be 
rectified  by  iteration.  Similar  expressions  are  defined  for 
6e(t)  and  6a(t).  The  time  series  of  radar  data  minus  the  reference 
trajectory  gives  the  three  data  residuals  6R(t),  6E(t),  and  6A(t). 
The  deviation  Ar(t),  defined  as  6r(t)  -  5R(t),  is  still  linear  in 
(X  ) ;  similarly  defined  are  Ae(t)  and  Aa(t).  Under  the  assumption 
of  a  normal  distribution  in  the  errors,  a  likelihood  expression 

L(X  )  is  then  defined  as  the  sum  of  the  squares  of  these  residuals 
n 

divided  by  their  respective  standard  deviations.  Referring  to 

N  data  points  and  to  standard  deviations  ct  .  a  ,  and  a  for  the 

Pea 
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radar  observables  we  have: 


t 

N 

L( X  )  =  -  Z  {[Ar(t)/a  ]2  +  [Ae(t)/o  l2  +  [ Aa( t) /a ]2 }* 

n  t=t  e  a 

o 

Maximization  of  L(X^)  yields  a  set  of  6  normal  equations  to 

determine  (X  ) . 

n 

The  procedure  used  in  the  computational  program  based  on 
the  analysis  is  as  follows: 

1.  A  data  set  of  N  points  is  given  along  the  trajectory. 

A  smaller  set  m  of  these  points  is  fitted  in  r,  e,  and  a  with 
polynomials  to  obtain  starting  initial  conditions. 

2.  The  equations  of  motion  are  integrated  with  these  ini¬ 
tial  conditions  to  obtain  a  reference  trajectory. 

3.  The  six  initial  conditions  are  then  perturbed  one  at  a 
time  and  new  trajectories  are  obtained.  Differences  in  r(t), 
e(t),  and  a(t)  are  formed  by  subtracting  the  reference  trajectory 
values  from  those  of  the  perturbed  trajectories.  These  differences 
are  fitted  with  polynomials  to  form  a  polynomial  coefficient  mat¬ 
rix  P(i,j,n).  The  subscript  i  from  1  to  3  refers  to  range,  ele¬ 
vation,  or  azimuth,  j  refers  to  the  degree  of  the  polynomial  up 


*The  inclusion  of  doppler  data  deviations  is  straightforward  but 
is  not  retained  here  for  simplicity. 
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to  a  maximum  of  J,  and  n  to  the  species  of  the  perturbed  initial 
condition.  Fourth  order  polynomials  are  used  for  range  residuals 
and  second  order  polynomials  for  the  angles. 

4.  Data  residuals  are  obtained  by  subtracting  the  refer¬ 
ence  trajectory  from  the  data.  These  data  residuals  are  fitted 
with  polynomials  to  form  the  coefficient  matrix  Q(i,j).  The  like¬ 
lihood  expression  L(Xn)  is  then  formed  as  follows: 


L(X  ) 
n 


1 


5.  The  partial  derivatives  of  L(Xn)  with  respect  to  the 

6  X's  are  set  to  zero  to  give  the  system  of  six  linear  equations 
for  determination  of  the  X's  which  in  turn  give  the  corrected 
initial  conditions  and  an  improved  trajectory.  If  the  resulting 
data  residuals  are  still  large  the  improved  trajectory  becomes 
the  reference  for  the  next  iteration. 

6.  As  the  iterations  continue  only  the  data  residual  poly¬ 
nomials  used  to  form  the  new  Q(i,j)  are  computed  while  the  recom¬ 
putation  of  P(i,j,n)  is  usually  not  necessary.  The  iterations 
cease  when  the  corrections  converge  to  negligible  values.  The 
RMS  residuals  of  data  minus  the  improved  trajectory  also  decrease 
with  every  iteration.  These  RMS  residuals,  rather  than  the  small 
corrections,  give  a  figure  of  merit  of  the  present  method. 
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In  the  different  stages  of  analysis  and  numerical  computa¬ 
tion,  different  coordinate  systems  are  adopted  either  for  logical 
reasons  or  for  convenience.  The  particular  coordinate  system 
used  in  the  analysis  of  the  likelihood  expression  and  initial 
condition  increments  consists  of  the  usual  range,  elevation,  and 
azimuth  for  position.  However,  for  the  velocity  part  of  the 
6-vector  giving  the  corrections  to  the  initial  conditions,  it 
has  been  found  useful  and  instructive  to  adopt  a  rectangular 
coordinate  system  of  special  orientation.  The  three  axes  point 
in  the  direction  of  (a)  the  range  rate  velocity,  (b)  the  remain¬ 
ing  component  of  the  velocity  perpendicular  to  the  range  rate, 
and  (c)  the  third  component  orthogonal  to  (a)  and  (b).  In  other 
words  the  Pg°  component  of  the  parameter  vector  is  defined  to  be 
zero.  This  choice  of  velocity  directions  is  useful  in  that  the 
range  rate  is  usually  known  quite  accurately  from  range  data  even 
if  doppler  data  is  not  available.  Thus,  the  velocity  part  of  the 
analysis  is  basically  a  determination  of  the  errors  in  the  two 
components  perpendicular  to  the  range  rate. 

III.  APPLICATION  TO  DATA 

To  illustrate  the  foregoing  procedure,  some  examples  of 
fitted  satellite  data  from  the  TRADEX  radar  will  be  given.  The 
computational  program  developed  uses  a  high-order  Runge-Kutta 
integration  and  a  gravitational  model  incorporating  spherical 
harmonic  terms. 
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The  range  and  elevation  radar  data  from  a  typical  satellite 
track  of  450  seconds  are  given  in  Figures  1  and  2.  The  data  were 
corrected  for  refraction  before  analysis.  After  three  iterations 
the  trajectory  residuals  in  range,  doppler*,  elevation  and  azi¬ 
muth  are  shown  in  Figures  3  to  6.  These  figures  illustrate  the 
quantitative  nature  with  which  the  derived  trajectory  follows  the 
data  up  to  quantities  all  of  the  order  of  the  noise  variations  of 
the  radar  observables.  This  is  particularly  remarkable  in  the 
case  of  the  range  where  the  departure  of  the  trajectory  from  the 
data  remains  within  the  order  of  25  feet  for  a  long  time  span  of 
400-500  seconds,  during  which  time  the  target  has  traveled  a  dis¬ 
tance  of  over  10  million  feet  and  the  range  varied  by  several 
million  feet.  This  case  is  typical  of  the  many  cases  analyzed, 

10  of  which  are  summarized  in  Table  1,  each  exhibiting  small  RMS 
residuals . 

A  study  was  made  to  determine  how  constant  the  trajectory 
initial  conditions  remain  as  shorter  data  spans  are  used.  An 
example  was  chosen  which  appeared  to  be  an  average  case.  First, 
290  seconds  of  data  (1  point  per  second)  were  used  and  the  initial 
conditions  obtained  for  this  trajectory  are  given  in  Table  2. 

The  velocity  components  given  are  a  standard  rectangular  set. 

Table  3  contains  the  velocity  differences  for  3  other  cases  using 

♦The  UHF  analog  doppler  data  is  just  shown  for  illustration;  it 
was  not  used  in  the  likelihood  expression.  The  sporadic  noise 
component  was  due  to  equipment  problems  and  has  been  eliminated. 
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shorter  data  spans.  These  velocity  differences,  which  can  be 
interpreted  as  errors  in  the  initial  velocity  estimates,  increase 
monotonically  as  the  data  span  is  shortened,  becoming  worst  for 
the  50-second  case.  On  the  other  hand,  Table  4  indicates  that 
the  values  for  the  RMS  residuals  remain  nearly  constant  for  all 
the  cases.  This  apparent  contradiction  is  explained  by  the  fact 
that  to  obtain  the  same  RMS  values  the  shorter  span  of  data  does 
not  require  as  much  accuracy  in  the  determination  of  the  initial 
velocity;  while  the  smaller  velocity  differences  obtained  for  the 
longer  spans  give  a  measure  of  the  accuracy  needed  in  the  velocity 
determination  in  those  cases.  Note  that  the  initial  range  rate 
receives  only  minor  changes.  Thus,  the  main  changes  are  in  the 
components  orthogonal  to  the  range  rate. 

We  shall  use  the  iteration  history  of  the  changes  in  the 
initial  velocity  for  two  cases  to  illustrate  some  remarks  about 
the  speed  of  convergence  and  the  special  velocity  coordinate 
system  mentioned  previously.  Table  5  contains  the  corrections 
in  the  initial  conditions  from  one  iteration  to  the  next  for  the 
first  case.  The  procedure  basically  converged  in  one  trial,  and 
this  fact  indicates  that  the  assumption  of  linearity  of  the  per¬ 
turbations  relative  to  the  corrections  is  a  valid  one.  The  case 
in  Table  5  is  one  for  which  a  large  component  of  the  initial 
velocity  is  in  the  range  rate  direction,  and  this  tends  to  sim¬ 
plify  the  determination  of  the  initial  velocity.  In  Table  6,  a 
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case  is  shown  for  which  the  initial  range  rate  is  a  very  small 
component  of  the  initial  velocity.  In  this  case,  the  procedure 
needed  three  iterations  for  convergence  but  again  the  largest 
corrections  were  in  the  velocity  components  orthogonal  to  the 
range  rate  velocity.  For  further  illustration,  the  second  and 
third  iteration  trajectories  of  this  last  case  will  be  given  to 
show  how  sensitive  the  range  fit  is  to  velocity  errors.  Figure  7 
contains  the  range  residuals  after  the  second  iteration  and 
Figure  8,  the  residuals  after  the  third  iteration.  The  addi¬ 
tional  velocity  correction  introduced  after  the  second  iteration 
is  only  about  1.5  feet  per  second;  yet  this  change  accounts  for 
the  difference  in  the  two  fits.  The  RMS  residual  values  for  the 
two  cases  are  given  in  Table  7.  The  determination  of  the  initial 
velocity  to  an  accuracy  of  one  foot  per  second  implies  a  very 
high  accuracy  in  determining  the  initial  angle  rates.  Since 
angle  rate  is  not  a  measured  quantity,  ordinary  angle  smoothing 
methods  could  not  derive  these  rates  with  an  accuracy  comparable 
to  the  global  fitting  technique  used  here. 

IV.  REFRACTION  CORRECTION  ANALYSIS 

Consideration  was  given  to  ascertaining  the  importance  of 
applying  refraction  corrections  to  elevation  data,  especially  at 
low  angles.  An  attempt  was  made  to  determine  if  the  refraction 
corrections  could  be  obtained  as  residuals  from  an  analysis  of 
uncorrected  data.  The  results  of  this  effort  were  quite 
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satisfactory.  Shown  in  Figure  9  are  the  elevation  residuals 
obtained  from  uncorrected  data  in  a  case  which  contained  low 
elevation  angles.  In  the  same  figure,  tabulated  refraction 
corrections  are  plotted,  and  these  values  essentially  agree  with 
the  residual  curves.  Additional  results  are  shown  in  Figure  10 
from  another  satellite  track.  For  the  cases  analyzed  which 
began  at  low  elevation,  it  was  found  necessary  to  start  the 
integration  from  a  high  elevation  to  obtain  satisfactory  conver¬ 
gence.  This  fact  coupled  with  the  residuals  obtained  when  no 
refraction  corrections  are  applied  point  out  the  necessity  for 
applying  refraction  corrections  to  radar  data  to  obtain  accurate 
trajectories . 

A  set  of  closed  formulas  for  refraction  correction  was 
derived  for  the  elevation  data.  The  essentials  of  the  deriva¬ 
tion  are  contained  in  the  appendix.  The  numerical  results  agree 
with  the  tabulated  values. 

V.  CONCLUSIONS  AND  DISCUSSION 

1.  A  technique  to  obtain  accurate  maximum  likelihood  tra¬ 
jectories  has  been  developed  and  applied  to  radar  satellite  data 
with  excellent  results. 

2.  The  measure  of  the  goodness  of  the  trajectory’s  fit  to 
the  data  lies  on  one  hand  in  the  final  RMS  values  of  the  residuals 
in  range  and  angles  and  on  the  other  hand  in  the  requirement  that 
the  fit  prevails  over  as  long  a  time  span  as  practicable.  The 
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typical  horizon-to-horizon  track  is  of  the  order  of  500  seconds, 
and  for  these  spans  RMS  values  of  25  feet  in  range  and  0.05 
degrees  in  angle  have  been  obtained. 

3.  The  best  initial  conditions,  especially  in  velocity, 
are  obtained  after  a  number  of  iterations  to  determine  the  cor¬ 
rections  to  these  conditions.  In  practice,  one  to  three  itera¬ 
tions  were  found  to  be  necessary  for  convergence,  and  this  de¬ 
pended  on  whether  the  target  was  initially  approaching  head-on 
(faster  convergence)  or  in  a  broadside  orientation  (slower). 

4.  The  repeated  iterations  essentially  obtain  corrections 
to  the  two  velocity  components  perpendicular  to  the  range  rate 
velocity.  The  magnitudes  of  these  components  can  be  of  the  order 
of  hundreds  of  feet/second  in  the  first  iteration  down  to  one 
foot/second  on  the  final  iteration. 

5.  It  is  essential  to  apply  the  refraction  corrections  to 
elevation  angle  data  in  cases  of  low  elevation  angles.  In  addi¬ 
tion,  if  the  target  appears  first  at  low  elevation,  it  is  found 
necessary  to  start  the  integration  from  the  high  elevation  angle 
end  to  obtain  reasonable  convergence. 

This  study  and  the  availability  of  the  related  computational 
program  presents  a  number  of  possible  extensions.  One  area  of 
study  would  be  to  extend  the  trajectory  to  the  next  pass  of  the 
satellite  and  compare  it  with  actual  radar  observations.  Such  a 
comparison  would  be  extremely  informative  and  is  currently  being 
developed . 
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The  permissible  approximation  of  the  18  influence  functions 
of  time  by  a  finite  polynomial  coefficient  matrix  P(i,j,n)  with 
less  than  200  elements  suggests  that  further  study  should  be  made 
to  determine  if  these  P(i,j,n)  elements  can  in  turn  be  approxi¬ 
mated  by  simple  functions  of  the  initial  conditions.  This  in 
turn  would  suggest  the  exploration  of  applying  the  method  to 
real-time  trajectory  determination. 

Another  extension  could  be  to  study  the  determination  of  a 
ballistic  coefficient  curve  which  is  described  by  a  few  parameters 
just  as  this  method  has  been  applied  to  finding  the  6  initial 
conditions  describing  a  target's  trajectory. 
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APPENDI X 


APPROXIMATE  FORMULAS  FOR  REFRACTION  CORRECTIONS  IN  AN  EXPONENTIAL 

ATMOSPHERE 

Refraction  corrections  in  elevation  under  a  single  exponen¬ 
tial  refractivity  model  are  usually  obtained  by  careful  numerical 
integration  and  tabulated  as  a  function  of  range  and  elevation. 

It  is  fairly  cumbersome  to  incorporate  this  table  into  a  computer 
program.  We  therefore  have  derived  a  set  of  approximate  formulas 
for  the  elevation  correction  as  a  function  of  range,  elevation 
and  the  two  parameters  a  and  B  defining  the  exponential  refrac¬ 
tivity  model. 

1.  Notation 

R  =  Range,  in  units  of  earth's  radius 
e  =  Apparent  elevation 
Y  =  True  elevation 

H  =  Altitude,  in  units  of  earth's  radius,  as  a  function  of 
R  and  Y 

n  =  Radio  refractive  index  =  1  +  a  exp  (-BH) 

Ae  =  e-Y  =  Refraction  correction  for  elevation 

2.  Formula  I:  For  elevation  <  4.5° 

i  2 

Ac(radians)  =  a  (Cos  Y(1  +  Sin  Y/R)(ttB/2)2  exp(B  Sin  ilr/2) 

| erf [  (B/2)  2(R  +  Sin  i[r)]  -  erf[(B/2)^  Sin 
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3. 


4. 


5. 


Formula  II:  For  elevation  >  4.5° 

Ae(radians)  =  a(Cos  ^  +  Cos  it  Sin  ifr/R) 

{[1-(1/B)(1  +  1/B)  /Sin2  it  ]  /S  i  n  ty 
-  exp  ( -BH)  [  l-(  1/B)  ( 1  +  1/B  +  H)/(R  +  Sin  it)2]/ 
(R  +  Sin  Y)j  +  a[exp  (-BH)  -1]  Cos  Y/R 
To  obtain  a  smooth  function  of  the  elevation  that  will  bridge 
over  the  transition  at  4.5°,  one  could  use  a  linear  interpola¬ 
tion  formula  between  I  and  II  for  the  interval  of  say, 

4°  <  e  <  5°. 

The  above  set  of  formulas  has  an  error  of  less  than  one  per¬ 
cent  as  compared  to  accurate  numerical  integration  results. 
Since  the  formulas  have  been  derived  in  terms  of  the  unknown 
true  elevation  Y,  an  iterative  procedure  has  to  be  used, 
starting  with  e  in  the  place  of  Y  in  the  formulas.  Two 
iterations  are  found  to  be  sufficient. 
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SUMMARY  OF  FITTED  TRAJECTORIES 


RMS 

C/3 

d) 

LO 

CD 

l> 

l> 

rH 

CO 

00 

s: 

03 

CD 

CO 

o 

O 

C2 

O 

CO 

CO 

CO 

4-5 

in 

C- 

rH 

rH 

in 

in 

O 

in 

!> 

00 

fcJ3 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

s 

03 

*H 

N 

c 

TJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

72 

s 

Otf 

c 

o 

•H 

4-5 

rt 

> 

0) 

rH 

w 


OOG^OOCOCOrHOCO  rH 


OCDCOOOrHCOOOO 
(NrHCNrHCOmCOrHCO  CO 


72 

s 

cr; 

0-5 

03 

03 

b£ 

03 

C 

a 

cn 

CO 

CO 

CO 

Tj< 

o 

Tf 

CO 

CD 

03 

O 

CO 

CO 

rH 

CO 

o 

rH 

o> 

in 

CO 

03 

rH 

rH 

CO 

CD 

rH 

CO 

00 

u 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

fcJD 

03 

TJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a  TJ 


72 

03 

s 


co 


o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

in 

C2 

CO 

CD 

CO 

m 

in 

in 

in 

in 

CO 

Tt< 

CO 

CO 

4-3 

CO 

03 

H 


C  CQ 

C0^mcDI>00C>OO  O 

cococococococo^o  o 

CD  CD 


15 


TABLE  II 


c o 

a 

!5 

o 

o 

w 


o 

LO 

CD 

N  0) 

• 

>  \ 

LO 

+-> 

o 

<H 

00 

rH 

o 

o 

LO 

CD 

• 

>>  CO 

LO 

>  \ 

+-> 

<H 

00 

CO 

CO 

O 

00 

o 

CD 

• 

X  (0 

CO 

CM 

>  \ 

o 

+-> 

rH 

O 

CH 

CO 

55 

rH 

CO 

CO 

Eh 

CO 

w 

Eh 

cc; 

o 

CO 

55 

O 


Q 

55 

O 

o 


c 

o 

•H 

+■> 

C3 

> 

O 

rH 

w 


CD 

CD 

Jh 

bi| 

CD 

T3 


o 

CO 

00 

m 

CM 


X! 

CD 

LO 

4-> 

CD 

U 

G) 

s 

bo 

• 

•rH 

CD 

o> 

N 

CD 

< 

CM 

•J 

h 

rH 

< 

CD 

o 

• 

HH 

rH 

CD 

CD 

Eh 

O. 

CO 

CD 

HH 

a 

\ 

CO 

J5 

o 

+-> 

00 

HH 

Q 

<H 

1 

CD 

O 

b£5 

+-> 

C 

CD 

o 

CCS 

CD 

PS 

rH 

LO 

16 


TABLE  III 


o 

55 


a 

a) 

l> 

CO 

o 

N 

CO 

CM 

o 

o 

> 

• 

• 

• 

<3 

>p 

O 

CO 

18 

a 

0) 

m 

rr 

>> 

CO 

00 

Oi 

o 

> 

\ 

• 

<3 

+-> 

o 

rH 

CO 

1 

rH 

co 

CQ 

Q 

W 

55 

hH 

2 

PQ 

O 

CO 

55 

O 


Q 

55 


CO 

55 

<3 

ft 

CO 

< 

H 


a 

o 

CM 

in 

CM 

X 

CO 

O 

a> 

> 

\ 

• 

<3 

+-> 

rH 

CM 

t> 

I 

i 

rH 

0) 

CM 

0) 

rH 

CO 

ft 

Sh 

o 

o 

o 

<3 

bo 

o 

o 

o 

0) 

• 

• 

• 

T3 

o 

o 

o 

O  Q 

55 

0) 

a> 

CM 

a 

O 

CD 

CO 

00 

C D 

ts 

hH 

N 

Sh 

o 

o 

O 

J  w 

<3 

bp 

o 

o 

O 

<3  H 

II 

<3 

CD 

• 

• 

• 

M  f£ 

TJ 

o 

o 

O 

H  O 

a 

1 

1 

4- 

hH  W 

hH 

55  CO 

<3 

hH 

55 

a 

M 

a) 

rH 

O 

•tS 

CO 

• 

• 

• 

CO 

<3 

\ 

o 

rH 

o 

w 

+-> 

1 

o 

<H 

55 

W 

35 

w 

ft 

ft 

05 

-p| 

CO 

00 

CO 

HH 

<3 

«h| 

CM 

1 

ft 

t 

17 


RMS  VALUES  FOR  CASES  CONTAINED  IN  TABLES  II  AND  III 


c 

0 

•H 

CD 

+-> 

CD 

03 

> 

b D 

CD 

0 

rH 

T3 

w 

CO 

CM 

O 

O 


LO  ^ 

t>  CO 

CO  ^ 

o  o 

o  o 


in 

c d 

o 

o 


0 

3 

0) 

e 

u 

•rH 

bD 

N 

CD 

< 

T3 

cd 

co 

t> 

o 

o 


CM  CO 

rH  00 

I> 

o  o 

o  o 


CD 

in 

o 

o 

o 


0) 

b£ 

C 

03 

es 


0) 

(D 

'-H 


o 

cd 


<T>  CO 

m  co 

rH  CM 


CM 

CD 


55 1 


in 

C 

o 

o 

CD 

c/) 


o 

G) 

CM 


o  o 

o  o 

CM  rH 


O 

m 


18 


TABLE 


W 

O 

3 


+-> 

0 

.Q 

a 

<N 

oo 

rH 

o 

o 

a) 

rH 

O 

o 

o 

o 

w 

• 

• 

• 

• 

• 

\ 

o 

o 

o 

o 

H 

+-> 

<N 

1 

> 

<3 


I* 

Eh 

hh 

O 


• « 


o 

a 

in 

o 

o 

00 

0 

0) 

00 

rH 

o 

o 

o 

w 

p> 

W 

• 

• 

• 

• 

l> 

> 

— 1 

\ 

+-> 

in 

<N 

o 

o 

o 

rH 

> 

i 

H  W 
M  CO 
15  <d 

HH  O 

55  W 

M  Eh 

CO  P§ 

w 


►« 

05 


15  O 
<  55 

a 

0) 

rH 

CO 

o 

o 

«  < 

•03 

V) 

rH 

rH 

o 

o 

in 

O  « 

> 

\ 

• 

• 

• 

• 

CD 

<1 

+-> 

o 

o 

o 

CO 

O 

1 

rH 

co 

h-i 

in 

15 

o 


w 

Eh 


C3 

O 

•H 

-P 

Ctf 

P 

0) 

+-> 

hH 


<N 


CO 


>» 

p> 

•H 

a 

o 

rH 

Q) 

> 


a 

p> 

o 

Eh 


19 


TABLE  VI 


Ol 

o 

H 

co 


both 

o 

in 

l> 

O 

CM 

o 

(D 

CM 

CM 

O 

o 

o 

03 

• 

• 

• 

• 

• 

p> 

\ 

CO 

O 

rP 

o 

o 

-1 

-P 

4H 

to 

rP 

CM 

<S 

01 

g 

>< 

H 

w 

O 

O 

*J 

W 

> 

-3 

i— i 

Eh 

HH 

55 

M 

55 

t— I 


w 

CO 

<d 

o 

w 

Eh 


•01 

o 

rH 

rP 

o 

CM 

o 

0 

CD 

O 

o 

o 

o 

-p 

m 

• 

• 

• 

• 

• 

H 

\ 

CM 

CO 

o 

o 

CO 

■p 

03 

1 

1 

1 

CM 

> 

<*J 

4H 

O 

t> 

CO 

CM 


CO 

w  w 
o  o 


55  55 

o 

m 

pP 

00 

o 

o 

•01 

CD 

m 

i> 

rP 

CO 

o 

o 

O 

> 

\ 

03 

CM 

o 

o 

o 

< 

-p 

CM 

rP 

1 

CO 

Pp 

1 

in 

w 

55 

o 

M 

Eh 

2 

W 

Eh 

HH 


c 

o 

•H 

^P 

rt 

P 

Q) 

P> 

w 


O  rP 


CM  CO 


p> 

•P 

o 

o 

rP 

CD 

> 


c3 

P> 

0 

Eh 


20 


c 

0 

0 

•H 

o 

4-> 

u 

a 

bJD 

> 

0 

0) 

T3 

r— t 

w 

00  *> 

CO  CO 

t>  t> 

o  o 

o  o 


XI 

XI 

xs 

0 

CO 

t> 

4-> 

0 

m 

LO 

3 

Sh 

i> 

l> 

C/0 

s 

bfl 

o 

o 

•H 

0 

• 

• 

PS 

N 

73 

o 

o 

2 

C 

CO 

^  w 

O  CO 
w  < 

M 

Eh  O 

hH 

<! 

> 

w  w 

W  Eh 

w 

Eh  <rj 

XI 

HH  P2 

CQ 

< 

CQ  W 

H 

o 

P5  Jz; 

o  < 

0 

rH 

!> 

Ph  W 

bto 

4-> 

CD 

CO 

0 

• 

• 

CO 

0 

o 

w 
>— 1 

pc: 

CO 

CM 

t— ' 

XI 

«3 

> 

CO 

s 

PS 

21 


ACKNOWLEDGMENT 


The  authors  wish  to  acknowledge  the  contributions  of  others 
to  this  work.  K.  Ralston  was  one  of  the  designers  of  the  TRADEX 
real-time  tracking  system  and  contributed  helpful  discussions  and 
computer  routines  to  this  work.  H.  M.  Jones  and  J.  P.  Penhune 
also  contributed  helpful  discussions.  V.  L.  Lynn  and  B.  Nanni 
were  helpful  in  supplying  the  satellite  track  data  from  the  TRADEX 
radar . 


22 


REFERENCES 


1.  J.  R.  Bauer,  W.  C.  Mason  and  F.  A.  Wilson,  "Radio 

Refraction  in  a  Cool  Exponential  Atmosphere",  TR-186, 
Lincoln  Laboratory,  (27  August  1958). 

2.  H.  Schneider,  "Multidimensional  Parameter  Estimation  by 
the  Summed  Weighted  Least  Squares  Minimization  of 
Remainders,"  with  Addendum  by  H.  Schneider  and  A.  R. 

Warwas,  Journal  of  the  Astronautical  Sciences,  (Fall  1964). 


3.  F.  C.  Schweppe,  "Algorithms  for  Estimating  a  Re-entry 
Body's  Position,  Velocity  and  Ballistic  Coefficient  in 
Real  Time  or  for  Post  Flight  Analysis,"  1964-64,  Lincoln 
Laboratory,  M.I.T.,  (3  December  1964). 


4. 


I.  I.  Shapiro,  "Prediction  of  Ballistic  Missile  Trajec¬ 
tories  From  Radar  Observations,"  McGraw  Hill  (1957). 


23 


RANGE  (10  ft) 


Figure  1:  Range  versus  time  for  test  33. 
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Elevation  angle  residuals  versus  time  for  test  33. 
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Azimuth  angle  residuals  versus  time  for  test  33. 
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Figure  7:  Range  residuals  versus  time  before  convergence  for 
test  600A . 
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Figure  8:  Range  residuals  versus  time  at  convergence  for 
test  600A. 
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ELEVATION  DEVIATION  (deg) 


Figure  9:  Elevation  angle  residuals  as  compared  with  tabulated 
refraction  corrections  for  test  40. 
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ELEVATION  DEVIATION  (deg) 


Figure  10:  Elevation  angle  residuals  as  compared  with  tabulated 
refraction  values  for  test  600B. 
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